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Anharmonic  lifetime  of  H and  D on  the  Si  surface 

I.  P.  Ipatova,  O.  P.  Chikalova-Luzina,  and  K.  Hessf 
Ioffe  Physico-Technical  Institute,  St  Petersburg,  Russia 
t Beckman  Institute,  The  University  of  Illinois,  USA 

Abstract.  The  study  of  H and  D adatoms  on  the  surface  of  Si  by  STM  technique  and  by  optical 
measurements  indicate  that  anharmonic  interaction  of  adatom  localized  vibrations  with  surface 
phonons  of  Si  defines  both  the  absorption  line  width  and  the  yield  of  the  desorption.  The  anharmonic 
lifetime  of  H and  D adatoms  is  calculated. 


STM  studies  of  the  Si  surface  terminated  by  H or  D demonstrated  that  the  lifetime  of  D- 
adatom  is  much  shorter  then  the  lifetime  of  H-adatom.  At  the  low  level  of  the  excitation,  the 
desorption  yield  demonstrates  the  temperature  dependence  which  indicates  on  the  essential 
role  of  the  H or  D vibrational  anharmonic  lifetime  in  the  process  of  desorption  [1]. 

The  study  of  1R  absorption  by  stretching  localized  vibrations  of  H adatoms  on  the  surface 
Si  (111)  [2]  has  demonstrated  the  anomalous  temperature  dependence  of  the  integrated 
absorption.  The  line  intensity  (integrated  area)  decreases  by  20%  as  the  temperature  is 
raised  above  the  temperature  130  K to  560  K.  Authors  of  the  paper  |2]  believe  that  the 
anomalous  temperature  dependence  follows  from  the  strong  anharmonic  interaction  of 
H-Si  stretching  modes  with  H-Si  bending  mode. 

All  the  facts  of  the  case  have  something  in  common  with  those  which  were  widely 
discussed  in  60-es  in  connection  with  optical  properties  of  U-center  localized  vibrations 
in  alcali  halides.  It  had  been  shown  by  Hughes  [3]  and  Ipatova  et  al.  [4]  that  anomalous 
temperature  dependence  of  the  integrated  intensity  can  be  understood  to  be  due  to  the  large 
difference  of  the  localized  frequency  and  the  frequencies  of  the  bulk  modes  of  alcali  halides 
(adiabatic  approximation).  It  had  been  shown  that  the  strength  of  the  absorption  by  the 
localized  mode  has  the  exponential  temperature  dependence  of  the  Debye- Waller  type. 

This  paper  deals  with  the  theory  of  anharmonic  line  width  and  integrated  absorption 
by  localized  vibrations  of  H-  or  D-adatoms  on  the  surface  of  Si  (1 1 1).  The  anharmonic 
interaction  of  Si-H  stretching  localized  mode  with  the  substrate  vibrational  Si  modes  is 
taken  into  account.  Since  the  Si-H  stretching  mode  frequency  is  <wioc(H)  — 2085  cm-1, 
stretching  mode  frequency  of  Si-D  is  1 5 1 6 cm- 1 [5],  and  the  maximum  phonon  frequency 
of  Si  is  col  (Si)  = 5 14  cm-1 , the  adiabatic  parameter. 


^loc 

holds  for  both  H-Si  and  D-Si  adatoms.  It  means  that  due  to  anharmonic  interaction  the 
slow  motion  of  Si-atoms  occur  in  an  averaged  field  of  fast  vibrations  of  light  adatom.  It 
is  shown  that  there  appear  static  displacements  of  Si  atom  equilibrium  sites.  Vibrations  of 
the  Si  atoms  occur  about  new  equilibrium  positions. 

The  shift  of  equilibrium  positions  results  in  the  exponential  temperature  dependent 
factor,  Debye-Waller  factor,  in  both  the  integrated  absorption  and  in  the  spectral  line  width 
of  the  light  absorption  by  localized  vibrations  of  the  H or  D adatoms.  This  exponential 
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factor  depends  on  the  temperature  and  on  the  mass  of  adatom.  It  is  different  for  the  localized 
vibrations  of  isotope  complexes  Si-H  or  Si-D. 

We  assume  that  the  adatom  (H  or  D)  of  the  mass  M'  is  situated  on  the  (1 1 1)  surface  of 
Si.  Adatom  is  supposed  to  be  bound  by  the  force  constant  in  the  “on-site”  configuration  to 
the  surface  atom  of  Si  at  the  lattice  site  l — (7,  ly  lz)  = (0  0 0)  = 0.  The  position  of  adatom 
is  / = (0  0 1 ) = 1 . The  mass  of  the  Si  atom  is  M. 

The  coefficient  of  the  absorption  of  1R  light  at  frequencies  co  ~ a>ioc  has  the  form  [6] 

e2  p Ea  Eg 

K (co)  = -cd  — Im  D^p  (1,  1;  co)  ^JL , (2) 

Cd  being  the  surface  concentration  of  adatoms,  i>o  is  the  volume  of  the  crystal  primitive 
cell  and  DRo  (1 , 1 ; to)  is  the  one-particle  retarded  Green’s  function  of  the  adatom  localized 
vibrations. 

The  retarded  Green's  function  ( 1 , I ; oj)  can  be  calculated  by  the  temperature  Green’s 

function  technique.  'I'he  temperature  Green’s  function  Dap(  1 . 1 ; ico„)  with  co„  =4 jtnT /h 
satisfies  Dyson  equation 

Dap(\,\-Uon)  = Dfp(l,\\iaon) 

+ Y!  iMn)  n yS(ls,  l' s':  ico„)  D$p(l's',  1;  icon),  (3) 

Is, l' s' 


where  Is,  l's'  are  positions  of  the  crystal  lattice  sites;  ny^(/.v.  Vs':  icon)  is  the  polarization 
operator  produced  by  anharmonic  terms  in  the  potential  energy  of  the  system.  Taking 
into  account  the  square  symmetry  of  the  adatom  position  on  the  surface,  one  can  find 
Daft  (1 , 1 * i (On  ) = D ( 1 , 1 ; i (On )&afi  and  v\ap  ( 1 , 1,  i COn  ) = n ( 1 , 1 , i(On)$afi • 

Adiabatic  parameter,  Eq.  (1),  allows  to  neglect  in  the  sum  over  (Is,  Vs')  in  Eq.  (3) 
non  diagonal  terms  £>(1,  Is;  ico„)  and  to  find  D ( 1 . 1;  icon)  from  algebraic  equation.  The 
retarded  Green’s  function  DR(  1,  Ecu)  is  the  analytic  continuation  in  a > of  temperature 
Green’s  function  D(\,  1 :i(o„).  One  can  separate  in  11(1,  1 ; co)  the  real  and  the  imaginary 
parts,  n(l,  l;o>)  = A(l,  I;  co)  — /T(l,  1;  co),  which  are  in  charge  of  the  anharmonic 
shift  Ren  = A ( 1 , I ; a>)  of  localized  frequency  and  the  decay  of  the  localized  vibration 
Imn  = r(\,  l;o>). 

The  quantity  F ( 1 , I ; oj)  is  defined  by  the  processes  which  satisfy  the  conservation  laws 
of  the  energy  and  the  momentum.  When  the  third  order  anharmonic  term  is  taken  into 
account,  in  the  first  approximation  of  the  perturbation  expansion,  F(l,  1;  cu]0C)  does  not 
vanish  only  when  col  < (Oioc  < 2 col-  The  condition  holds  neither  for  H-  nor  for  D- 
adatoms.  In  order  to  find  T(l,  1 ; co)  the  anharmonic  terms  of  the  4th  order  should  also 
be  taken  into  account.  There  are  no  limitations  of  this  type  for  the  shift  A(l,  1;  oj):  all 
possible  virtual  anharmonic  phonon  processes  contribute  to  A.  As  a result,  A(o>)  > F (oj) 
and  the  shift  A is  defined  by  the  third  order  anharmonisity.  The  expansion  of  11(1,  1 ; oj) 
near  co  « cuioc  gives 

H 

M7 

where 


Dk(\,  Ecu)  = 


A/(<W]oc) 


)=T 


coc(T)  + iy(co)' 
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Here  AVioc)  = (%A) 


V dco2  ) 


loc 

Calculations  of  polarization  operator  n(l,  1;  to)  and  corresponding  coc  and  y with  the 
help  of  the  temperature  Green’s  function  technique  allow  to  use  the  parameter  Eq.  (1)  and 
to  find  the  main  contribution  to  K (co) 


K (co)  — t'4 


hv o \2M 


t rujoc  / 


-2  W(T)_ 


-2 W(T)  y 


(o>  — ay.)2  + (e  2W(r)r)2 


Here  the  quantity  W(T)  is  defined  by  the  3rd  order  anharmonic  terms 


2 W(T)  = Y] — — ( ) 

\2M'co\oc) 


*(3)  2 CNx+J). 
{Tuox)2 


(6) 


(7) 


where  cox,  N,_  are  the  normal  mode  frequencies  and  the  occupation  phonon  numbers  of 
vibrations  in  semi-infinite  Si,  <t>®  is  the  closure  of  the  anharmonic  force  constants  with 
the  phonon  polarization  vectors. 

It  follows  from  Eq.  (6)  that  the  quantity  e~2wri  1 r(&j|oc)  is  the  line  width  of  optical 
absorption  by  high  frequency  localized  mode.  At  low  temperatures  T < ©d,  ©d  being 
Debye  temperature,  the  factor  e~2w  does  not  depend  on  temperature:  e 2"/(0)  = const.  At 
high  temperatures  T > ©d  the  Debye- Waller  factor  is  the  linear  function  of  temperature. 

Two  contributions  to  the  quantity  T are  known:  the  relaxation,  V j , and  the  dephasing, 
T2-  At  low  temperatures  T < ©d,  the  main  contribution  to  the  line  width  for  H adatom 
T(H)  is  shown  to  come  from  dephasing  mechanism  of  broadening.  1^2  (H)  is  vanishing  as 
> at  j o.  In  case  of  D-adatom,  the  dephasing  T i defines  the  line  width  in  the 
temperature  interval  T > T\  only.  At  T = 7\  ~ 60  K.  the  relaxation  broadening  T i (D) 
becomes  more  important. 

At  low  temperatures  T < ©d,  the  comparison  of  the  line  width  for  H-  and  D-adatoms 
shows  that  at  T > Ti  30  K the  main  mechanism  of  broadening  is  dephasing  F2  for  both 
H-  and  D-adatoms.  Order  of  magnitude  evaluation  shows  that  F2(H)  > F2(D).  At  T < 73 
the  relaxation  broadening  Fi(D)  becomes  more  important  for  D-line  width  and  H-line 
width  becomes  less  than  D-line  width:  IMH)  < F i (D).  At  high  temperatures  T > ©d, 
the  temperature  dependence  of  the  broadening  is  defined  by  T(T)  e~2  W(T>  where  T(T) 
is  dephasing  process  for  both  H-  and  D-adatoms.  Numerical  evaluations  shows  that  at 
T > 250  K,  the  relaxation  between  the  line  widths,  r(D)e_2w(D)  > r(H)e-2w(H), 
holds  due  to  exponential  factor. 

The  integrated  absorption  is  equal  to 


, e 

J — I dco  K (<w)  — q 

Jo 


n 

2 M'  ft)]0c 


e-2  W(T)_ 


(8) 


It  follows  from  Eq.  (8)  that  the  integrated  intensity  J has  exponential  temperature 
dependence  which  is  stronger  in  case  of  H-adatom  than  for  D-adatom.  Order  of  magnitude 
evaluation  has  shown  that  /(H)  2.7(D)  what  is  in  the  good  agreement  with  experimental 

data  from  [6]  where  /(H)  = 1.9.7(D)  is  obtained. 
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